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Abstract 

We investigate the magnetic properties of small transition metal clusters us- 
ing a simple statistical model, which requires some input data from ab initio 
spin-density functional calculations. In our study, we consider a thermo- 
dynamically equilibrated ensemble of clusters with different structures, spin 
multiplicities, and ground state energies. We calculate the physical proper- 
ties of this system by weighting the individual configurations according to 
the Boltzmann statistics. We find that presence of isomers with very similar 
ground state energies, yet very different magnetic properties, gives rise to a 
rich magnetic behaviour of the system which differs significantly from what 
would be expected for single configurations. We apply the present model 
to determine the magnetic susceptibility of a cluster ensemble of Langevin 
paramagnets. Our results show that some of the anomalies in the magnetic 
behaviour of transition metal clusters might be understood in the framework 
of our model which is, of course, limited by the extremely high computational 
effort needed to obtain the input data. 
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I. INTRODUCTION 



The magnetic properties of small transition metal clusters have been of growing interest 
in the past few years. However, there has not been any study about the dependence on 
temperature so far. 

As shown in a large variety of papers (e.g. Refs. [0,0) small clusters exhibit, unlike 
the bulk, smooth structural transformations, which one might call isomer hopping, and 
which occur between pure solid and fluid phases over a relatively wide temperature range, 
in a so-called coexistence phase. In the following we will be concerned with the effect of 
those properties on the magnetic behaviour of clusters. The basic idea is that a cluster 
of a specific size might have two or more structures with different magnetic moments and 
that these structures or isomer states occur with their statistical probability.^ As shown 
below, this simple assumption causes a strong deviation of the paramagnetic behaviour of 
the magnetic susceptibility from the Curie law. As an interesting feature, a dependence 
on the strength of an external magnetic field might occur, even at relatively small fields. 
Transition metal clusters seem to be natural candidates for the occurrence of such effects 
because of their different possible electron configurations with very similar energies, which 
in turn yield quite different average magnetic moments. 

II. THEORY 

We consider to have n isomers with magnetic moments /ij and ground state energies Eo(i), 
calculated using an appropriate spin-density functional Hamiltonian. In the presence of an 
external magnetic field the energy changes. First order quantum mechanical perturbation 
theory yields 



In the following we will view different configurations, even if they have the same symmetry and 
only different lattice constants, as different isomers. 
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= g^^^MB,^, + ^B^, + vD ' (1) 

where g is the Lande factor and the magnetic field is pointing in the 2;-direction (see Refs. 
IP^). The first term describes a paramagnet, which will be considered exclusively in what 
follows and can be described with a purely classical description p, as 

£;^ = -^-H. (2) 

Now we are able to write down the partition function as 

n I 

Z = V / dcos9exp{-pEo{i) + pfiiHcos(9)) 
i -^-1 

n n 

= exp{-pEo{z))—— smh{Pfi,H) . (3) 

Indeed this is an extremely simple partition function, which can be improved to any com- 
plexity by adding all degrees of freedom of the clusters, e.g. by substituting the summation 
over the isomers by integration over the whole configuration space over the atomic positions 
or by adding a summation over different spin multiplicities of the clusters. From Eq. 
the average magnetic moment (/x) and the magnetic susceptibility x can be calculated from 
Eq. (D easily [|], using 

_ Eiexp(-/5Eo(i)) cosh(/3/iiiJ) 1 



E,/ir'exp(-/3Eo(0)sinh(/?/i,if) /?H ' 

X=limi^, (5) 

H^Q VH ^ ' 

if Nmjy is taken as the number of particles per unit volume.0 In the case of small magnetic 
fields /i H ^ /ebT we get 



^ The thermal property related to x by the dissipation-fluctuation theorem is 7 = — ■^^hhIh.Z' 
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n 1 

(/.) ^ Zo E exp(-/3Eo(^))- /3 /x,^ (6) 

i 

and 

/V " 1 

X ^ -^Zo Y.^M-f^E,{i))- p /i? . (7) 

Here, Zq{[3) is the partition function defined in Eq. (|^) for zero magnetic field. Eq. 
reveals simply that x is a linear superposition of the susceptibilities of the individual isomers 
weighted by their thermal probability. 



III. RESULTS AND DISCUSSION 

The general behaviour of the paramagnetic susceptibility according to Eq. (|^) is plotted 
in Fig. 1 for a hypothetical system with only two major isomer states with magnetic moments 
/ii and /i2 and their energy difference AE = Eo(2) — Eo(l). Unhke a purely paramagnetic 
behaviour according to the Curie law, the plots reveal local minima and maxima. A similar 
behaviour has recently been experimentally found by Cowen et al. for Fe28 clusters in 
supercages of NaY Zeolite. The locations of the extrema only depend on the ratio between 
the magnetic moments /ii/yU2 and AE. For a given ratio ^i/ the position of the minimum 
depends in an almost linear fashion on AE, as can be seen in Fig. 2. The extraordinary 
sensitivity of these results on AE is illustrated by the fact that a change of AE by only a 
few meV moves the minimum by hundreds of degrees Kelvin. 

In Fig. 3 the average magnetic moments of the canonical ensemble are plotted for various 
magnetic fields. These plots reveal the fact, that in the presence of a magnetic field the 
occupation probability of the isomers changes dramatically due to the additional magnetic 
energy. 

If one views H to be small in the sense of Eq. (5), the magnetic susceptibility can be 
easily calculated. For varying magnetic fields this calculation yields a behavior which is very 
similar to that of a varying ground state energy difference (see Fig. 4). 
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At first sight, configurations of different symmetry - such as bcc versus fee struetures - 
of Fe clusters seem to be the best candidates to show the effects discussed above. Even if the 
energy difference between two isomers might be very small, the transition from one isomer 
to the other may involve a complex concerted motion of atoms and may be associated with 
a nonzero activation energy. Thus a thermal relaxation might take quite a long time, and a 
hysteresis might be observable. 

Lee and Callaway have found interesting results for Crg and Vg bcc clusters, which 
are given in Table I. With varying lattice spacing the average magnetic moment changes by 
a factor of up to five, whereas the ground state energy changes are only ^ 0.04 Ry/atom. 
This is especially interesting since even a simple spatial expansion of the transition metal 
clusters intuitively is more probable than a structural transformation. 

From our results we infer that it is extremely hard to obtain phenomenological results for 
the magnetic behaviour of clusters from spin-density functional methods, since the results 
depend very sensitively on ground state energy differences, which are often very small and 
comparable in magnitude to the precision of these ab initio methods. In addition, for larger 
clusters, the number of relevant isomers increases dramatically. 

Nevertheless, we have shown that in the model case of two cluster isomers, the rich 
behaviour of the magnetic susceptibility as a function of temperature and the external 
magnetic field is a sensitive tool to probe structural properties of clusters. On the other 
hand, once the general magnetic response due to the above discussed effects is understood, 
one might use the susceptibility of the clusters as a very sensitive thermometer. 

Moreover, by applying a sufficiently strong magnetic field, one might be able to suppress 
one of the transition metal cluster isomers, provided the clusters have ample time for struc- 
tural rearrangement. An extended study of these effects, including specific examples, will 
be published elsewhere. 
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TABLES 

TABLE L Average magnetic moments /x (/XB/atom) and total energy per atom Etotai (Ry) for 
Crg and Vg clusters, as listed in Ref. 0. 
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FIGURES 

FIG. 1. Magnetic susceptibility X / A*B ^ for an ensemble consisting of two magnetic config- 
uration states {ill =1.0 //B , M2 =10.0 hb)- AE is the energy difference between the two states. 

FIG. 2. Location of the local minima of the magnetic susceptibility as a function of AE for 
various ratios /X1//X2 of the magnetic moments. 

FIG. 3. Average magnetic moments at various magnetic fields H, for AE = 0.01 eV, m = 
1.0 hb, and /X2 = 10.0 /xb- 

FIG. 4. Magnetic susceptibility %/ hb ^ at non-zero magnetic fields H (eV/ /xb), for AE = 
0.01 eV, Hi = 1.0 hb and 112 = 10.0 hb- 
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